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Use of Chebychev Polynomials in Thin Film Computations 

Klaus D. Mielenz 

(August 17, 1959) 

From Herpin's expression for the mth power of a multilayer matrix, very simple closed 
formulas are derived for the matrices and optical constants of any multilayer with a periodic 
structure. 

According to Epstein's theorem, any symmetrical multilayer is equivalent to a fictitious 
monolayer. A simple expression for the equivalent index and thickness of this monolayer is 
deduced for the case of a periodic and symmetrical sequence of equally thick films. 

As compared to any other method of numerical computation, the suggested formulation 
provides a considerable saving of time and work. In a numerical example, this saving 
amounts to about 80 percent. 



1. Fundamental Relationships [1, 2] l 

The electromagnetic field at the plane of entry to a 
multilayer as in figure la is determined by 
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are the complex index and the optical thickness in 
phase units of the v\\\ Layer (n P = refractive index, 
'■c v = absorption coefficient, fx v = permeability, d v = 

hysical thickness, X= wavelength, and <fc,=angle 

f incidence). 
The amplitude transmission and reflection coeffi- 
cients of the multilayer are 
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with 



C=M E +H , D=M E -H . 
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Various methods for computing R and T have 
been suggested in the literature. They all have the 
disadvantage of being based upon recurrence rela- 
tions that make it necessary to calculate the desired 
quantities in a cumbersome stepwise manner. Al- 
though such a procedure seems to be inevitable in 
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Figure 1: Optical multilayers. 

(a) General case and denotations, (b) periodic, and (c) periodic-symmetrical 
multilayer. (In figures H> and c, the Individual layers a and b do not necessarily 

represent single films.) 



general, a much simpler approach is possible with 
"periodic" and "periodic-symmetrical" multilayers 
as represented in figures 1 1> and c. Since such multi- 
layers are of considerable importance in thin film 
work, the computation method developed hereafter 
is of great practical significance. 

2. "Periodic" Multilayers 

A multilayer in which the same sequence of films 
is repeated twice or more often is a "periodic" 
multilayer. According to Herpin's theorem [A] any 
multilayer, and therefore the fundamental period of 
layers as well, may be expressed as a fictitious bilayer 
the matrix of which we shall call 9f a 2l & . If the 
period occurs m times the matrix eq (1) reads 
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3. 'Teriodic-Symmetriccd" Multilayers 

A multilayer for which the indices and thicknesses 
are the same as encountered from either side, 

is a " symmetrical" multilayer. It can be replaced 
by a fictitious monolayer (Epstein [4]). 

If such a multilayer consists of [m + (l/2)] times 
a fundamental period, 



(t) =(w ^(A> 

it is a "periodic-symmetrical" multilayer. 



(6) 



4. The mth Power of a Matrix of Unity 
Determinant 

A simple closed expression for (2l a 2t & ) m is the key 
to eqs (5) and (6). Herpin [5] has shown that the 
powers of a four-element matrix can be expressed 
by Lucas polynomials, and Abeles [6,1] has observed 
that these are reduced to Chebychev polynomials 
if the basic matrix is of unity determinant (which 
is the case for multilayer matrices) . This principle, 
however, has not been developed further since. 
Let % be any four element matrix of determinant 
unity and write 

\a 21 a, 2 / 

+d<r 3 , a 2 -(b*+c*+d*) = l, (7) 

with 
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(Pauli spin matrices) . Then, we obtain 

2P=2a2I-ao. (9) 

Following Herpin, we set 

W»=S m -MS m - 2 <r , (10) 

whence we arrive at 
?t 



Putting 



- 1 =SJL-S m . 1 c =fSr , n=S M . 1 fP-S m - a i 
=2aS m - 1 %-S m - 1 <r -S m - 2 %. 

2a=aii-\-a 2i =X, 



(11) 



we find by comparison of coefficients, as recurrence 
formula for the SJs, 



S m (X) — XS m _i (X) — S m _2 (X) . 
Equation (9) yields the initial values, 
S {X) = 1, S 1 (X)=X, 
and then (12) leads to 



sjao-'S (-D-C") 
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(12) 
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where [m/2] denotes the largest integer contained 
in m/2, e.g., [5/2] = 2. Explicitly, we have 



$o=l, 
S 1 =X, 

S 2 =X 2 -1, 
S 3 =X 3 -2X, 
Si=X<-3X 2 +l, 
S 5 =X 5 -4X 3 +3X, 
S 6 =X 6 -5X 4 +6.Y 2 -1, 
S 7 =X'-6X 5 +lOX z -iX, 
S s =X s -7X li +15X i -lOX 2 +l, 
etc. 



(14a) 



The S m 's defined by these equations are the 
Chebychev polynomials of the second kind [7], 



S n (X) = 



sin (ra + l)9 



or 
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X=2 cos 0, (15a) 
A"=2 cosh $. (15b) 



Another form is 
1 



s m (xy- 
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-(X-TpP=l) m + 1 ]. (16) 



For real arguments, X=x, these polynomials are 
also real (even though this is not obvious in eq (16) 
for |#|<2). If, in eqs (15), and $ are to be real 
for reasons of convenience, (15a) must be used for 
|x|<2, and (15b) for |z|>2. 

The desired matrix (2l a 2l&) m ma y n( >w be obtained 
as follows: Form 

^*=(yT) and X=a n +a 22 . (17a) 

\tt 2 i 0-22/ 

Then, find S m -i(X) and S m - 2 (X) and write 

(nr=( Sm -% n ~ Sm - 2 Q Sm ~ ia i \ d7b) 
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5. Application to Multilayers 

According; to (5), (6), and (17), we have 

F P =S m - 1 (X).F ab -S m -*(X).F„ (18) 



and 



Fps — Sm-1 (X) 'Faba — Sm-2(X) -F a , 



(19) 



where F stands for either field strength, E or H, and 
because 4 of (4) for the quantities C and D as well. 
The subscripts P and PS refer to the whole periodic 
and periodic-symmetrical multilayers, respectively; 
the subscripts s, a, ab, and aba refer to the uncoated 
substrate 4 , the bottom mono-, bi-, and trilayers, 
respectively. 

Thus, any of these quantities can be expressed as 
a simple linear combination of the corresponding 
quantities of much simpler multilayers. 

By forming the matrix product in (16), one finds 
the argument of the Chebychev polynomials in (18) 
and (19), 

Mj+Mj . n . D ,„. 



X=2 cos B a cos B b 



6. Layers of Equal Thickness 

The equations of section 5 are still further simpli- 
fied if each film in the multilayer has the same 
optical thickness, 

(21) 



B a= B t -H. 



Then, we have 



MaM* 



(22 



While for the periodic multilayer the mathematical 
formulation itself is not simplified, a very simple 
formulation is obtained for the periodic-symmetrical 
multilayer: 

With (21) and (22), the matrix of the trilayer (aba) 
can be shown to be 



njm a =x^ a -^u\ 



where 



Hi" l =( 
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(23) 



(24) 



is the inverse of ?l„. Then, (17) and (12) show that 



m a nt,) m ^a=s m (X)% l -s m - 1 (X)^ 1 - 



(25) 



hat 



Because of (13), (23) is contained in (25) as the par- 
ticular case m=l. 

Thus we see that, if all films are equally thick, it 
is not necessary to express the periodic-symmetrical 
multilayer in terms of the bottom trilayer and the I is a real number. 



uncoated substrate, as in (19). Instead, it can very 
simply be expressed in terms of the two basic matrices 
% a and 2l & , without any need for multiplication of 
mat rices whatsoever. 

Epstein's theorem [4], according to which any 
symmetrica] multilayer is equivalent to a fictitious 
monolayer, was already mentioned. The index M m 
and the thickness B m of the monolayer corresponding 
to a periodic-symmetrical multilayer with equally 
thick films may now be obtained as follows: 

Consider (25) and write 



(8a3W w 2l«= 



(cos B n 
%M m si 



Mj m M (2G) 



%M m sin B m cos H„ 



Then, comparison of coefficients yields, for the 
thickness B m and the index M m of the fictitious 
monolayer, 



cos B m =[S m (X)-S m . 1 (X)] cos B, 



M*= 



M a SJX)+M b S m - 1 (X) 



M a Sm{X)+ M b Sm - l(X) 



(27) 



(28) 



(According to Epstein [4], the signs of B m and M m 
have to be chosen such that B TO ->0 for \-^oo, and that 
always Re (M m )>0.) 

7. Dielectric Multilayers 

The above formulation constitutes a considerable 
simplification of practical computations. 

An automatic computer can provide for itself 
the needed Chebychev polynomials by computing 
them according to eqs (14) or (16), regardless of 
whether AT is a real or a complex number. 

For desk calculations, however, numerical values 
of Chebychev polynomials of complex arguments 
cannot be found except with rather complicated 
calculations. This leads us to the restriction that 
X always should be a real number which will be 
true only if the multilayers are purely dielectric, and 
if the individual matrices 2t a and 2t & represent individ- 
ual films. (The replacement of multilayers, even if 
purely dielectric, by fictitious mono- or bilayers may 
yield complex indices.) These assumptions lead to 
the important class of alternating dielectric layers, 
for which we have 



M a =m a 



M b =m b =n b , 



(29) 



and (for equally thick layers thai are quarter wave 
films at a wavelength X ) 



p 2X 



X=x=2 



B a =B b 

(n tt +nt,) 2 



n a n b 



sin 2 /3 



(30) 
(31) 
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For such real arguments, the Chebychev poly- 
nomials may be found with the aid of: 
J (a) Numerical tables: 12-decimal values of the 
first 12 £ m 's, for 0<x<2 with intervals 0.001 in x, 
have been published by the National Bureau of 
Standards [7] . Jones and co-authors [8] have published 
similar but less complete tables. For negative argu- 
ments, find S m (\x\) and use the relation 



S m (-x) = (-l)»S n (x). 



(32) 



(b) Equation (15) and tables of trigonometric or 
hyperbolic junctions: If, for instance, S 4 (3.745) is 
looked for, one finds from the table that 3.745 = 2 
cosh 1.24. Hence $=1.24, 5<£> = 6.20, and 

o /Q7i ,x sinh 6.20 246.37 

S4(3J45)= sinTri^ = T^83n =155 - 62 - 

(c) Direct computation: The S m 's may also be 
computed from eq (14) or (16). For small values 
of m it is easier to use eqs (14), but for large m's 
(which, however, rarely occur in multilayer work) 
eq (16) is faster. 

8. A Practical Example 

Consider a high reflection multilayer consisting 
of 1 1 alternating zinc sulfide and magnesium fluoride 
films on glass, 

n =l (air), n a =2.3 (ZnS), 
7i, 6 =1.38 (MgF 2 ), n s =1.52 (glass). 

Let all films be a quarter wave thick at X =5460.74 A, 
and compute the amplitude reflection coefficient R 
for X=4358.35 A, i.e., for 



7T 5460.74 



2 4358.35 



= 112° 45.86', 



sin 0=0.92211, cos £=—0.38694. 
Then we have according to (1), (24), and (25), 



f E o\ To / x/-0.38694 0.40092f\ 

UH^H 2.12085^ 



S 4 (x) 



0.38694 J 

0.38694 -0.66820i 
27252; -0.38694 
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From (31), we obtain x= — 1.62789, so that we may 
look up in reference [7]: S 5 (x) = 0.94004, S±(x) = 
0.07257. 
Thence, 



©=[("?: 



-0.36374 
99363; 



0.37688;\ 
-0.36374 / 



+ 



\0. 



02808 
09235; 



0.04849; 
0.02808 



)](!.,) 



■( 



-0.33566 
2.08603; 



0.42537;\/l ^ 
-0.33566 Al.52 y 



< 



-0.33566+0.64656;\ 
-0.51020+2.08603;/ 



Finally, (3) and (4) yield 



R-- 



0.17454 -1.43917; 
-0.84586 +2.73259;' 



Starting with the given n's and with (3, this result 
was obtained with 17 individual steps of multiplica- 
tion or division and 9 steps of addition or subtraction. 
Four numerical values had to be looked up in tables. 

In comparison hereto, it takes 96 multiplications, 
and 48 additions or subtractions, with 1 numerical 
value to be looked up, to arrive at the same result 
by means of the widely used recurrence method for 
admittances [9], This may be estimated to be about 
five times as much work. 
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Selected Abstracts 

Measurement of the aging of rubber vulcanizates, 

J. Mandel, F. L. Roth, M. N. Steel, and R. D. 
Stiehler, J. Research NBS 63C No. 2, Ul (1959). 

A study of aging data in the literature and of measurements 
made at the National Bureau of Standards indicates that 
ultimate elongation is the best of the tensile properties for 
characterizing the deterioration of rubber vulcanizates during 
storage at various temperatures. Ultimate elongation de- 
creases during aging for all types of rubber vulcanizates; 
whereas, tensile strength and modulus may increase, de- 
crease, or remain essentially unchanged. The change in 
ultimate elongation over prolonged periods of storage cannot 
be expressed by a simple mathematical equation. How- 
ever, during most of the useful storage life of a rubber vul- 
canizate, the elongation decreases approximately linearly 
with the square root of time. The data indicate that for 
some vulcanizates an estimate of storage life at room tem- 
perature can be made from measurements of ultimate elon- 
gation at two or more elevated temperatures. 

Excitation mechanisms of the oxygen 5577 emission 
in the upper atmosphere, E. Tandberg-Hanssen and 

F. E. Roach, J. Research NBS 63D No. 3, 319 (1959). 

Possible excitation median isms for the green 5577 emission 
are considered in the light of recent data on the dynamics 
of the upper atmosphere. Photochemical reactions as affect- 
ed by mass motions as well as excitation directly due to the 
mass motions are analyzed. 

Method for measuring local electron density from 
an artifical satellite L. R. O. Storey, J. Research 
NBS 63D No. 3, 325 (1959). 

A method is proposed for measuring the electron density at 
known points in the outer ionosphere, by the use of vlf 
receiving equipment in an artificial satellite, in conjunction 
with a vlf transmit fcer on the ground. The transmitter would 
radiate continuous waves, which would be propagated 
through the ionosphere in the 'whistler' mode. The basis 
of the method is a measuremenl of the local wave admittance 
of the medium, by comparison of the signals received on an 
electric dipole and on a loop. 

A further proposal is made for an integrated vlf satellite 
experiment, in which several different types of observation 
would be made simultaneously. 

Reflectors for a microwave Fabry-Perot interfer- 
ometer W. Culshaw, IRE Trans, on Microwave 
Theory and Techniques, Vol. MTT-7, No. 2, 221 
(1959). 

The advantages of microwave interferometers for wave- 
length and other measurements at millimeter wavelengths are 
indicated, and a microwave Fabry-Perot interferometer dis- 
cussed in detail. Analogous to the cavity resonator, this 
requires reflectors of high reflectivity, small absorption, and 
adequate size. Stacked dielectric plates, and stacked planar or 
rod gratings are shown to be suitable forms of reflectors, and 
equations for the reflectivity, optimum spacing, and band- 
width of such structures are derived. A series of stacked 
metal plates with regularly spaced holes represents a good 
design of reflector for very small wavelengths. Fringes and 
wavelength measurements at 8-mm wavelength are given for 
one design of interferometer, these being accurate to 1 in 10 4 
without any diffraction correction. For larger apertures and 
reflectors in terms of the wavelength, errors due to diffraction 
will decrease. 



The nature, cause, and effect of the porosity in elec- 
trodeposits, III. Microscopic detection of porosity, 

F. Ogburu and D. W. Ernst, Plating 46, 831 (1959). 

A new technique called "parallel sectioning" was used for 
investigating the porosity of nickel electrodeposits. Some 
of the coating is polished oil* in a layer parallel to the basis 
metal. After taking a photomicrograph of i his area, another 
layer of metal is removed. This procedure is continued until 
the basis metal has been reached. It is possible to construct a 
three-dimensional model of the coating. Various pore types 
are listed and illustrative photomicrographs are given. 

Mechanism of contraction in the muscle fiber-ATP 
system, L. Mandelkern, A. S. Posner, A. F. Diorio, 
and K. Laki, Proc. Nat. Acad. Sci. U.S. 45, No. 6 y 

814 (1959). 

The axial contraction that develops when glycerol-treated 
rabbit psoas muscle was immersed in ATP solutions of vary- 
ing concentrations was measured. The length-concentration 
diagram has the characteristics typical of a cooperative phase 
i ransformal ion and the shrinkage can be attributed to melting 
of the initially axially oriented polypeptide chains. This 
conclusion was substantiated by wide angle x-ray diffraction 
studies of the native and shrunken fibers which show the 
structural changes expected on melting. Thus, the same 
principles that govern the contractile processes in a wide 
variety of fibrous macromolecular systems are also operative 
in the muscle- ATP system. 

Microwave spectrum of methyl germane, V. W. 

Laurie, J. Chem. Phys. 30, No. 5, 1210 (1959). 

The </ = 0— >1 and 1— >2 transitions of 28 isotopic species of 
CH 3 GeH 3 have been measured. From the rotational con- 
stants obtained, the following structural parameters have 
been calculated: r C H= 1. 083 ± 0.005 A°, r Ge u= 1.529 ±0.005A°, 
r CG e= 1.9453 ±0.0005 A°, <£HCH = 108°25' ±30', <£HGeH 
= 109°15 / ± 30'. The K=\ transitions of the asymmetrically 
deuterated species are split by internal rotation about the 
C — Ge bond. With the assumption of a threefold sinusoidal 
potential, an internal harrier of 1239±25 cal/mole has been 
determined from these splittings. Analysis of the Stark 
effect of several species gives a dipole moment of 0.635 ±0.00(> 
D. From observed hyperfme structure the nuclear (juad- 
rupole coupling constant of Ge 73 has been calculated to be 
+ 3 Mc. 



Measurement of ozone in terms of its optical ab- 
sorption, R. Stair, Advances in Chem,. Series of the 
Am. Chem. Soc. y No. 21, 269 (1959). 

The unique absorption spectrum of ozone provides an ideal 
physical basis for measuring its concentration in the atmos- 
phere even in the presence of significant quantities of other 
atmospheric pollutants, whether of gaseous or particulate 
character. Various types of optical equipment have been 
considered, both for measurement of the total amount of 
ozone and for determination of its vertical distribution and 
horizontal concentration. Natural sunlight furnishes a suit- 
able and convenient light source for measuring the total 
amount and vertical distribution of ozone. Special sources 
having high radiant intensity within the spectral region of 
2500 to 3600 A, where ozone has a high optical absorption, are 
desired for use in measuring the horizontal concentration of 
ozone The light source, whether the sun or some special 
source such as a mercury arc lamp, may be employed with a 
simple filter radiometer or with a more or less elaborate prism 
or grating spectroradiometer as desired. In most of the 
recent work at the National Bureau of Standards a double, 
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(jiiMi-i z prism spectroradiometer has been used at Washington, 
D.C., Climax, Colo, Los Angeles, Calif., and Sunspot, N. Mex., 
in ozone studies. 

Lower bounds for eigenvalues with application to 
the helium atom, N. W. Bazlev, Proc. Nat. Acad. 
Sci. U.S. 45, No. 6, 850 {1959)." ' 

Let ibea self -ad joint operator with domain D in a Hilbert 
space £). Suppose A = A' + A where A is self-adjoint and 
A' is positive definite. The eigenvalue problem for A, whose 
solution we assume known, gives rough lower bounds. If u% 
(*=1,_. . ., k) are k discrete eigenvectors of A and if P%= 
{A')~ l Ui (i—1, . . ., k) exist then one can substantially 
improve the lower bounds. The theory is applied to the 
helium atom operator. 

Spectroscopic evidence for triatomic nitrogen in 
solids at very low temperature, M. Peyron, E. M. 

Horl, H. W. Brown, and H. P. Broida, J. Chem. 
Phys. 30, No. 5, 1304 (1959). 

Additional wavelength, intensity and lifetime measurements 
of the radiation emitted from solid nitrogen containing 
trapped atoms have led to the interpretation that a weakly 
bound triatomic molecule, N2— N, is an emitting species. 
The three lowest electronic levels of atonic nitrogen 2 P, 2 D, 
and 4 S are involved in the eight line groups which have been 
found. Isotopic substitution has confirmed this model. 
Evidence also has been found for an N 2 — O molecule similar 
to the N 2 — N. 

Thermodynamic properties of helium at low tem- 
peratures and high pressures. D. B. Mann and 
R. B. Stewart, NBS Tech. Note 8 (PB151367) $1.25. 

This is a compilation and correlation of the present data on 
the thermodynamic properties of helium below 20° K. The 
existing (best) values are selected. The results are presented 
in the form of temperature-entropy and enthalpy-entropy 
diagrams. Pressures to 100 atm, temperatures from 0° K to 
20° K and specific volumes from 5 liters/kg to 800 liters/kg 
are presented. 

On the perturbation of the vibrational equilibrium 
distribution of reactant molecules by chemical re- 
actions, K. E. Shuler, 7thSymp. (Intern.) on Combus- 
tion, London and Oxford, Aug. 28 to Sept. 3, 1958, 
Combustion Inst. p. 87 (Butterworths Sci. Pub., 
London, England, 1958). 

A mathematical analysis is made of the stepwise excitation of 
an assembly of harmonic oscillators with an irreversible dis- 
sociation limit by a solution of the transport equations using 
Gottlieb polynomials. The variations in the distribution 
function as determined by values of the dissociation energy 
and the conditions for perturbation of the equilbrium Boltz- 
mann vibrational distribution because of chemical reactions 
are discussed. 

Other NBS Publications* 

Journal of Research, Section 63B, No. 2, October- 
December 1959. 75 cents. 

Applications of a theorem on partitioned matrices. 

E. V. Hay ns worth. 
Capacity requirement of a mail sorting device: 

II. A. J. Goldman. 
Analytic comparison of suggested configurations for 

automatic mail sorting equipment. B. K. Bender 

and A. J. Goldman. 



New method of solution for unretarded satellite 

orbits. John P. Venti. 
Effect of sudden water release on the reservoir free 

outflow hydrograph. Vujica M. Yevdjevich. 
Uniform asymptotic expansions for Weber parabolic 

cylinder functions of large orders. F. W. J. 

Olver. 

Journal of Research, Section 63C, No. 2, October- 
December 1959. 75 cents. 

Conical coaxial capacitors and their advantages. M. 

C. Selby. 
A photoelectric followup and recording system, and 

its application to remote observations of the beam 

in high precision balances. H. A. Bowman and L. 

B. Macurdy. 
A stroboscopic vibration analyzer. S. Edelman, R. 

Brooks, S. Saito, E. Jones, and E. R. Smith. 
Evaluation of lens distortion by the inverse nodal 

slide. Francis E. Washer and Walter R. Darling. 
Evaluation of lens distortion by the modified gonio- 

metric method. Francis E. Washer and Walter R. 

Darling. 
Proposed criteria for defining load failure of beams, 

floors, and roof constructions during fire tests. J. 

V. Ryan and A. F. Robertson. 
Conductive flooring for hospital operating rooms. 

Thomas H. Boone, Francis L. Hermach, Edgar H. 

MacArthur, and Rita C. McAuliff. 
Measurement of the aging of rubber vulcanizates. J. 

Mandel, F. L. Roth, M. N. Steel, and R. D. 

Stiehler. (See above abstracts.) 

Journal of Research, Section 63D, No. 3, Novem- 
ber — December 1959. 70 cents. 

Radio-refractive-index climate near the ground. B. 

R. Bean and J. D. Horn. 
Path antenna gain in an exponential atmosphere. 

W. J. Hartman and R. E. Wilkerson. 
Effect of atmospheric horizontal inhomogeneity upon 

ray tracing. B. R. Bean and B. A. Calioon. 
Correlation of solar noise fluctuations in harmoni- 
cally related bands. L. R. O. Storey. 
A monochromatic low-latitude aurora. F. E. 

Roach and E. Marovich. 
Pattern synthesis for slotted-cylinder antennas. 

James R. Wait and James Householder. 
Central Radio Propagation Laboratory exponential 

reference atmosphere. B. R. Bean and G. D. 

Thayer. 
Excitation mechanisms of the oxygen 5577 emission 

in the upper atmosphere. E. Tandberg-Hanssen 

and F. E. Roach. (See above abstracts). 
Method for measuring local electron density from 

an artificial satellite. L. R. O. Storey. (See above 

abstracts) . 



Thermodynamic properties of helium at low tem- 
peratures and high pressures, D. B. Mann and 
R. B. Stewart, NBS Tech. Note 8 (PB151367) 
$1.25 (see above abstracts) . 
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Transmission loss in radio propagation — II, K. A. 

Norton, NBS Tech. Note 12 (PB151371) $3.00. 
Technical considerations leading to an optimum 

allocation of radio frequencies in the band 25 to 

60 me, K. A. Norton, NBS Tech. Note 13 

(PB151372) $2.50. 
The role of Government research laboratories, A. 
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